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Abstract 

We estimate the decay rates of r/c 27, r/^ 27, and J ^ e'^e~, 
tp' — > e"'"e~, by taking into account both relativistic and QCD radiative correc- 
tions. The decay amphtudes are derived in the Bethe-Salpeter formaUsm. The 
Bethe-Salpeter equation with a QCD-inspired interquark potential are used to 
calculate the wave functions and decay widths for these cc states. We find 
that the relativistic correction to the ratio R = V{r]c — > 2'y)/T{J/ip — > e~^e~) is 
negative and tends to compensate the positive contribution from the QCD 
radiative correction. Our estimate gives T{r]c — >■ 27) = (6 — 7) keV and 
T{rj'^ — > 27) = 2 keV, which are smaller than their nonrelativistic values. The 
hadronic widths r(r/c ^ 2g) = (17 - 23) MeV and r(r?^ 2g) = (5 - 7) MeV 
are then indicated accordingly to the first order QCD radiative correction, if 
ctsim-c) = 0.26 — 0.29. The decay widths for bb states are also estimated. We 
show that when making the assmption that the quarks are on their mass shells 
our expressions for the decay widths will become identical with that in the 
NRQCD theory to the next to leading order of v'^ and a^. 



Charmonium physics is in the boundary domain between perturbative and non- 
perturbative QCD. Charmonium decays may pvovide useful information on under- 
standing the nature of interquark forces and decay mechanisms. Both QCD radiative 
corrections and relativistic corrections are important for charmonium decays, because 
for charmonium the strong coupling constant as{mc) ~ 0.3 [defined in the MS scheme 
(the modified minimal subtraction scheme)] and the velocity squared of the quark in 
the meson rest frame f ^ ~ 0.3, both are not small. Decay rates of heavy quarkonia in 
the nonrelativistic limit with QCD radiative corrections have been studied (see, e.g., 
refs. [1,2,3]). However, the decay rates of many processes are subject to substantial 
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relativictic corrections. In the present paper, we will investigate relativistic correc- 
tions to the pseudoscalar quarkonium decays such as r]c ~^ 27 and 27, and give 
an estimate of their widths by taking into account both relativistic and QCD radia- 
tive corrections. (For a brief report on this result, see also ref.[4].) For comparison 
we will also study the leptonic decays of the vector charmonium such as J/'^ — > e^e~ 
and ■0' — e+e~. 

These pseudoscalar charmonium decays are interesting. Experimentally, the branch- 
ing ratio of rj^ 27 may provide an independent determination of at the charm 
quark mass, but the measured T{r}c — >■ 27) ranges from 6 keV to 28 keV [5], and the 
measured i^c total width is also uncertain. As for the r]'^, its exsistence needs to be 
confirmed, and its two gamma decay mode is being searched for by the E835 exper- 
iment at the Fermi Lab pp collider, and its hadronic decay modes are being studied 
by BBS CoUaboration at BEPC. 

A lot of theoretical work have been done on charmonium and, in particular, on 
these pseudoscalar charmonium decays [4, 9- 14]. Nonrelativistic quark model gives 
r(r7c — 27) = %.bkeV (using the observed J/ij^ leptonic width as input), while the 
QCD sum rule approach predicts a value of 4.6 ± 0.4/ceK[12]. 

Recently, there have been significant progresses in the study of heavy quarkonium 
decays based on a more fundamental approach of the NRQCD (nonrelativistic QCD) 
effective theory[15,16]. The factorization theorem was further discussed, and some 
important issues (e.g., the infrared divergences in the P-wave state decay rates) were 
clarified in this study. The NRQCD theory combined with nonperturbative lattice 
simulations have achieved many interesting results on heavy quarkonium spectrum 
and decays[17,18,19]. 

In this paper, we will use the Bethe-Salpeter (BS) formalism[20] to derive the 
decay amplitudes and to calculate the decay widths of heavy quarkonium. The meson 
will be treated as a bound state consist of a pair of constituent quark and antiquark 
(i.e., higher Fock states such as \QQg > and \QQgg > are neglected, which may 
be justified to the first order relativistic corrections of S wave heavy quarkonium 
decays) and described by the BS wavefunction which satisfies the BS equation. A 
phenomenological QCD-inspired interquark potential will be used to solve for the 
wavefunctions and to calculate the decay widths. Both relativistic and QCD radiative 
corrections to next-to-leading order will be considered based on the factorization 
assumption for the long distance and short distance effects. 

We first consider the rjc 2j decay. This process proceeds via the cc annihilation. 
In the Bethe-Salpeter (BS) formalism the annilination matrix element can be written 
as follows 

(0 \QIQ\P)^I d'qTr [I{q, P)xp{q)] , (1) 

where | P) represents the heavy quarkonium state, P{q) is the total (relative) momen- 
tum of the QQ, Xp{q) is its four dimensional BS wave function, and where /(g, P) is 
the interaction vertex of the QQ with other fields (e.g., the photons or gluons ) which, 
in general, may also depend on the variable (the time-component of the relative 
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momentum). If I{q,P) is independent of (e.g., if quarks are on their mass-shells 
in the annihilation), this equation can be written as 

(0 \QIQ\P) = j d'qTr[I{q,P)^p{q)], (2) 

where 

$P(?) = / dq\p{q) (3) 

is the three dimensional BS wave function of the QQ meson. Note that in this ap- 
proximation the decay amplitude is greatly simplified and only the three dimensional 

BS wave fTinction is needed (but this does not necessarily require the interquark in- 
teraction to be instantaneous). In the BS formalism in the meson rest frame, where 
Pi — —p2 = q, P = (M, 0), and Pi{p2) is the quark(antiquark) momentum, M is the 
meson mass, we have 

$^"(?) = AU^)7°(l + 7°)^7°A'-(-^)/(^), (4) 

where $p (q), and $p (q) represent the three dimensional wave functions of 0^ and 
1~ mesons respectively, ^ = e^j^, is the polarization vector of 1~ meson, (p and / 
are scalar functions which can be obtained by solving the BS equation for 0" and 1" 
mesons, and A+(A_) are the positive (negative) energy projector operators 

Aliq) = A+(pi) = -L(£; + . + ^/), 

Al{-q) = A4p2) = ^(E - 7O7 . P2 - ^7°), 

E = yj^ + m?. (5) 

For process r]c — > 27 with the photon momenta and polarizations qi, ei and 52, £2, 
the decay amplitude can be written as 

T = (0 I cr^.(g)c I 77e)e^(Ai)e2lA2) + (0 | cV'^Mc I r/c)e^(A2)e^(Ai), (6) 

where Pi{p2) is the charm quark (antiquark) momentum, p = Pi — qi, p' = Pi — q2, tti 
and M represent the masses of c quark and rjc meson respectively, and where 



^fiAq) = iiJi4-^iu, r^^(g) = lui-^r^iu, (7) 

p — m ^ p' — m 
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for Q = c. Since p? +P2 = ^1 ^ usual we take[l,13] 



= (8) 



Thus, — \M — = 0, = \M — g° = 0, the amphtude T becomes independent 
of q^. Employing Eqs. (2) and (4), we get 

where 

B = B', 

B = 'J^J 7^2 X L2 ^2 ^-irV^W- (10) 



(^^ +m2^ (^g^ + ql +w? -2q-qi^ 

Using gi • ei = and q2 • €2 — 0, it is easy to get the decay width 

r(7;e^27) = 3MWe^|5|% (11) 
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In the nonrelativistic (NR) hmit (M ~ 2r?7, q ^ 0) 



where we have used the relation 



j ip{q)dq^^^{% (13) 

where -0(0) is the Schrodinger wave function at origin in coordinate space. Substitut- 
ing (12) into (11), we get 

r^^(r;e - 27) = I27:a'e% \m? (14) 

where T^^lrjc — * 27) represents the decay width of 77^ 27 in the nonrelativistic 
limit, which is consistent with that given in ref. [1]. The QCD radiative correction 
to this process has been given in ref.[l]. Recently, in the framework of NRQCD the 
factorization formulas for the long distance and short distance effects were found to 
involve a double expansion in the quark relative velocity v and in the QCD coupling 
constant q;s[15]. To the next to leading order in both v"^ and ag, as an approximation, 
we may write 

Tin. 27) = SMWeJ \Bf (1 - M^^), (15) 

TT 



where the strong coupling constant as{mc) is defined in the MS scheme (the modified 

-^2 

minimal subtraction scheme). By expanding B in (10) in terms of q /m , to the 
next to leading order of we have 



16i /■ 11 ? 

B 



M{M^ + 4m2) 
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(16) 



We see that the relativistic kinematic effect is to suppress the rjc — > 27 decay width. 

For comparison with the process J/ij) ^ e^e" , wc also give the decay amphtude for 
the QQ annihilate into an electron with momentum ki and helicity ri and a positron 
with momentum /c2 and helicity r2- Here the interaction vertex /(P, q) = — ie7^, 
which is independent of q'^, and the amplitude can be written as 

T = e'cQiO I I J/^)u^^{ki)rvM^. (17) 
Define the decay constant fv by 

fvMe, = (0 I C7^c \J/^)^Jdq Tr[^,^^(q)], (18) 
where e^^ is the polarization vector of J/ijj meson. Then with (4) we find 



where E = \l (I +m?. Summing over the polarizations of the final states and aver- 
aging over that of the initial states, it is easy to get the decay width 

V{J/^l^ e+e-) = ^7ra^e|/^/M. (20) 

In the nonrelativistic limit it is reduced to the well known result 

r^i?( j/^ ^ e+e-) = IGTra^el |V'(0)|' /M^. (21) 
Including also the QCD radiative correction [1], we will get 

r(W - e-^e-) = ^™^e^f (1 - M^). (22) 
To the next to leading order of v^, Jv is expressed as 

Again, the relativistic kinematic correction is to reduce the leptonic decay width. 
Comparing the two photon width with the leptonic width, we get 

_ r(77e ^27) _ 9 3 , \B\' ^. (^s{m,) 

and in the nonrelativistic limit it becomes 

^ = r^i^(J/V; ^ e+e-) " 3^^ + l-^^^^)- (25) 
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In fact, there are two sources of relativistic corrections : 1) the correction of rela- 
tivistic kinematics which appears exphcitly in the decay amphtudes; 2) the correction 
due to inter-quark dynamics (e.g. the weU known Breit-Fermi interactions), which 
mainly causes the correction to the bound state wave functions. In general, due to 
the attractive spin-spin force induced by one gluon exchange for the 0~ meson, the 
Tjc wave function at origin becomes larger than its nonrelativistic value, one might 
expect the width of rjc 27 to be enhanced after taking relativistic corrections into 
account. However, because the kinematic relativistic correction to the decay rates is 
in the opposite direction and can be even larger, the overall relativistic correction to 
the decay width of rjc is found to be negative. 

To calculate the decay widths, we need to know the wavefunctions (p{q) for the 
0~ meson and f{q} for the 1~ meson, which are determined mainly by the long 
distance interquark dynamics. In the absence of a deep understanding for quark 
confinement at present, we will follow a phenomenological approach by using QCD 
inspired interquark potentials including both spin-independent and spin-dependent 
potentials, which are supported by both lattice QCD calculations and heavy quark 
phenomenology, as the interaction kernel in the BS equation. We begin with the 
bound state BS equation [20] in momentum space 

iii - mi)xp(g)(^2 + ^2) = ^ / d^kG{P, q - k)xp{k), (26) 

where qi and q2 represent the momenta of quark and antiquark respectively, G{P, q—k) 
is the interaction kernel which dominates the interquark dynamics. In solving Eq.(26), 
we will employ the instantaneous approximation since for heavy quarks the interaction 
is dominated by instantaneous potentials. Meanwhile, we will neglect negative energy 
projectors in the quark propagators which are of even higher orders. We then get 
the reduced Salpeter equation[20] for the three dimensional BS wavefunction ^p{q) 
defined in (6) 

Mq) = po_i^_^/ +7° / d'kG{P, q- k)Mkh'^-: (27) 

where G{P,q~ k) represents the instantaneous potential. 

We employ the following interquark potrntials including a long-ranged confinement 
potential (Lorentz scalar) and a short-ranged one-gluon exchange potential (Lorentz 
vector) [21] 

V{r) = Vs{r) + j^®rVv{r), 
(1 - e~°"') 

ar 

4 0- ir) 

Vv{r) = --^e~-\ (28) 
6 r 
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where the introduction of the factor e""'' is to regulate the infrared divergence and 
also to incorporate the color screening effects of the dynamical light quark pairs on the 
QQ linear confinement potential[22]. In momentum space the potentials become[21] 

G'(P)=G'5(P)+7m®7^G'v(P), 
Gs{P) = --S'{P) + , 

= - J,^. (29) 

where as{P) is the quark- gluon running coupling constant and is assumed to become 
a constant of as P — > 

o^siP) = ^ i . (30) 

^' ln(a + P /AIcd) 

The constants A, a, a and Aqcd sue the parameters that characterize the potential. 

Substituting (4) and (29) into Eq.(27), one derives the equation for the 0~ meson 
wavef unction (p{q) in the meson rest frame 

= —^{-{Eg,Eg2 + m^m2 + f) [ d'KGsiq- k) - 4Gv{q- k))^i{k) 
-{Eg,m2 + Eg2m,) [ d^k{Gs{q- k) + 2Gv{q- fe)) !?'!!?' V^i(fc) 

HE,. + E..)jd^kGsiq- m- 

+ (mi - ma) / (fkiGsiq- k) + 2Gy{q- k)){q- fe) ^ 7 f'' ^S)^ (31) 



where Egi = ^J(p + mf, Eki = yk'^ + mf, {i = 1, 2), and 

{mi+m2 + Eqi + Eg2){Egim2 + Eg2mi) 
= 4fl..g^(m.+m.) 

The normalization condition / d^qTr{^^ {q)^{q)} = {2n)~^2M for the BS wavefunc- 
tion leads to [21] 
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For the 1 meson we have 

Mh{q) = {E,r + E,,)h{q) 



Egi + mi + Eq2 + m2 r , , -e^ 



J m-i + mo 



/■^3>/^ or< r^^^{Eki + mi){Ek2 + m2) + P 
/ d k{Gs{q -k)- 2Gv{q - k)) 

J + £/fc2 + mi + 7712 

-gfc2 mi + -gfci 77l2 

mi + m2 

+(Egi7772 + ^52^1) / d - k) — — — ■ ■ fi[k) 

J Ekl + Ek2 + 7771+7712 

-(mi+7772) / d^k{Gs{q-k)+AGv{q-k)){q-k)fi{k) 

-2{E,i - E,2) I d'k{Gs{q- k) - 2Gy{q- k)){q ■ k) J^^'l"^' h{k) 

J tiki + 7771 + -C'fe2 + "7.2 



+ f d'k{4Gs{q-k) -SGviq-km-kf-^ +/r , ^ 

J hiki + 7771 + ^k2 + ^2 

-2(7771 - 7772) / d^kGs{q- k f''~^^\ q- k)fi{k) 

J mi + 7r72 

+{E,i + 3^,2) / d'kGsiq- k){q- k)fi{k) 

-{<6E,i + 2£;,2) / d'kGviq- k)q- kfi{k)}, (34) 



where 

,^ Eqi+mi + Eg2 + m2 ,^ , . 

fM = — T^r-ET /(€)■ (35) 

The normahzation condition / d^qTr{(^\q)(^{q)] = {2n)~^2M for the BS wavefunc- 
tion leads to [21] 

To the leading order in the nonrelativistic limit, Eqs.(31) and (34) are just the or- 
dinary nonrelativistic Schrodinger equation with simply a spin-independent linear 
plus Coulomb potential. To the first order of v^, Eqs.(31) and (34) become the well 
known Breit equations for the 0~ and 1~ mesons with both spin- independent and 
spin-dependent potentials from vector (one-gluon) exchange and scalar (confinement) 
exchange. 

For the heavy quarkonium cc and bb systems, rrii = m2 = m, Eqs.(31) and (34) 
become much simpler. By solving these equations we can find the wave functions for 
the 0~ and 1~ mesons. Here not only the ground state (IS) wave functions but also 
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the first radial excitation wave functions (2S) are obtained. They are shown in Fig.l 
and Fig. 2. 

Substituting the obtained BS wave functions into (10), (15), and (19), (22), re- 
spectively, we then get the decay widths for both 0~ and 1~ charmonium states. In 
the calculation following parameters have been chosen 

mc = l.hGeV, A = 0.23^6^^ Kqcd = 0.18Ge\/, 

a = QmGeV, a = e = 2.7183. (37) 

With these parameters the 25" — 15" spacing and J/ip — rjc splitting are required to fit 
the data. We then get 

r(?7c ^ 27) = 6.2keV, r{r]'^ 27) = l.SkeV, 
TiJ/ijj e+e-) = 5MeV, r(V'' ^ e+e") = 2.7A;ey". (38) 

Our results are satisfactory, as compared with the Particle Data Group experimental 
values[5] T{t]c 27) = 7.0t?:?fceV, r(J/^ ^ e+e') = 5.36 ± 0.29A;eV, r(V'' ^ 
e+e~) = 2.14 ± 0.21A;ey. Here in above calculations the value of as{mc) in the QCD 
radiative correction factor in (15) and (22) is chosen to be 0.29 (refs.[3,14]), which is 
also consistent with our determination from the ratio of B[J/i/j 3g) to B^J/t/j — > 
e+e") (see refs.[4,9]). 

In order to see the sensitivity of the decay widths to the parameters especially the 
charm quark mass, we have also used other two sets of parameters 

mc = lAGeV, X = 0.2AGeV^] 

mc = l.GGeV, A = 0.22GeV^, (39) 

with other potential parameters {Aqcd, Oi, a) unchanged, and found 

r(r]c 27) = 7.0(5.5)A;ey", r(r]'^ 27) = 1.7(1.5)A;ey (40) 

for rric = 1.4 (1.6) GeV, where the experimental value of T{J/ip — e+e") = 5.36keV 
(as input) and the calculated ratio R in (24) are used to give predictions for the 
pseudoscalar decay widths. 

We see that for smaller charm quark masses T{ric 27) gets enhanced. This 
tendency is in line with the QCD sum rule result [12]. Our estimate that T{r]c ~^ 27) = 
{6-7)keV is consistent with the GLEG data[6] T{ric 27) = (5.9l?i±1.9)A;e1/, and 
the £;760 data[7] 7±3keV, and slightly smaller than the L3 data[8] (8.0±2.3±2.4)A;ey. 
Our results for rjc and 77^ distinguish them from the nonrelavistic values, which can 
be obtained by using the ratio R^^ (25) and the experimental values of r[J/t/j{t/j') — > 
e+e-] 

T^^ir], 27) = 8.5keV, T^^iv'c ^ 27) = SAkeV. (41) 

In particular, our prediction T{7]'^ — > 27) = 2keV is significantly smaller than its 
nonrelativistic value. 
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We may further use these results to give an estimate for the total widths of rjc and 
77^. Note the branching ratio 

BiP 27) ^ ^ = ^ X ( -^) (42) 

is free of the relativistic correction. Using the calculated two gamma widths r(P — > 
27) = &.2{l.%)keV for P — ric{ri'^ and the strong coupling constant at the mass of 
the charm quark as{m^ — 0.26 — 0.29[3,9,14], we will get 

rtot(^c) = 17-23 MeV, 

Ttotiic) = 5.0 - 6.7 MeV. (43) 

This is the prediction for the total widths up to the next to leading order of QCD 
radiative corrections, but higher order corrections may further modify this result. 
With the present Particle Data Group values[5] VtotiVc) — 10.3t|;4Mey and V{ric — > 
27) = 7 .Q'^^keV , however, a value of as{mc) ~ 0.20 will be indicated, which is 
significantly lower than expected from other experiments and theoretical studies on 
the QCD scale parameter. Therefore, it will be very interesting to see the accuracy of 
the experiment or to take the higher order QCD radiative correction more seriously. 

Moreover, for the hh states, with = A.^GeV and other potential parameters 
unchanged, we find 

r(776 ^ 27) = QAQ>keV, T{r]'^ 27) = Q.2lkeV, 
r(T ^ e+e-) = 1.36fcey", r(r ^ e+e-) = O.TSfceV. (44) 

Here ^^(mb) = 0.20(3,9,14] is used in the QCD radiative correctons. We see that the 
relativistic corrections become smaller for hh than for cc states. 

We finally discuss the relation between our approach and the NRQCD theory. In 
fact, our decay widths can be written in terms of the standard Schrodinger wavefunc- 
tion (with relativistic corrections) ipschio), which is related to (p{q) (or f{q)) through 
the normalization condition (33) (or (36)) which leads to 

{2nf J d'q rschiQ)^Sch{q) ^ 1. (45) 

For the above discussed pseudoscalar(P) and vector(V) heavy quarkonium decays (see 
(10), (15), and (19), (22)) to the next to leading order in v"^ and ag we then have 
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In above expressions M is the mass of the meson. In previous calculations we have 
taken M as their observed values for 1]^. and J/ip. However, we may take the on-shell 
condition, which assumes the quark and antiquark to be on the mass shell (see (8)) 



= = M/2 = E = ^Jm? + (48) 
to replace the observed value of the meson mass M then (46) and (47) will become 

r(P ^ 2,) ^ 1^(1 - /Ad - m 

r,, . . _ i^^^ip)),/ A,i - (50) 

It is easy to see that to the first order of v'^, in coordinate space (49) and (50) can be 
expressed as 



np - 27) = ^(1 - ^-^^^^mm' + ^,Re{R*m'Rm], (51) 



r{V e+e-) = ^(1 - i^^i|pl)[|i?(0)r + ^i?e(i?*(0)V^i?(0))], (52) 

where -R(O) is the Schrodinger radial wavefunction at the origin of the P [P = rjc, 7]'^) 
or V {V — J/ip,ip') meson. These expressions are exactly the same as that given in 
ref.[15] with the NRQCD effective theory, if we identify our bound state wavefunctions 
with their regularized operator matrix elements, i.e.: 



/27r 

= yYe-<0|xW|V^>, 



R{^)^\h^e-<Q\x^aiP\V>, (53) 



V^i?(0) = -^^e < 0|xV(^ T>fi^\V > [1 + Oiv^/c^)]. (54) 

In the NRQCD theory, the expectation values of the quark operators are well de- 
fined[15] and can be calculated with lattice simulations, which is a more fundamental 
method for describing nonperturbative dynamics than the quark potential model. In 
our approach the wavefunctions (and their derivatives) are estimated on the basis of 
the QCD-inspired potential model by solving the BS equation. Although this is not 
a first principle theory and it is difficult to control the systematic accuracy within 
the potential model, it may provide a rather useful estimate of the decay rates, since 
not only the zeroth order spin-independent potential but also the first order spin- 
dependent potential i.e. the Breit-Fermi Hamiltonian, which stems from one gluon 
exchange and has a good theoretical and phenomenological basis, are considered in 
the calculation, and different quark massess are also chosen to estimate the uncertain- 
ties in the calculation. In fact, the potentials are required to reproduce the observed 
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mass difference between r)c and J/ip and the J/ip leptonic decay width, and then give 
predictions for the pseudoscalar mesons. This may reduce the uncertainty in the cal- 
culation of pseudoscalar decay widths. Nevertheless, for more reliable estimates we 
hope that these decay widths of heavy quarkonium can be eventually calculated from 
more fundamental theoretical methods, e.g., the lattice QCD simulations. It will be 
interesting to see the numerical results in the NRQCD approach and compare them 
with our results. 

In summary, we have estimated the photonic widths and hadronic widths for 
pseudoscalar heavy quarkonium states, and the leptonic widths for vector heavy 
quarkonium states as well, by taking into account both relativistic and QCD ra- 
diative corrections. The photonic widths of r)c and r)'^ tend to take lower values than 
the nonrelativistic result. We hope that experiments especially the E835 and BES 
experiments or experiments at the Tau-Charm Factory in the future will be able to 
have more accurate measurements on the branching ratios and the total widths for 
the rjc and r]'^ particles. This will provide the basis for testing theoretical predictions. 
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very useful conversations on the BES and E760-E835 physics programs especially 
regarding the study for rjc and 77^ particles. This work was supported in part by the 
National Natural Science Foundation of China, and the State Education Commission 
of China. 
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